home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 090 / byt85jul.lbr / ORBIT.FQR / ORBIT.FOR
Text File  |  1985-09-15  |  4KB  |  91 lines
  1.  
  2.       SUBROUTINE ELEMNT 
  3. C    ROUTINE TO CALCULATE ORBITAL ELEMENTS OF THE MAJOR PLANETS
  4. C    FOR ECLIPTIC AND EQUINOX OF 1950.0
  5.       DOUBLE PRECISION OELMT(8,8),CLAT(8),CRAD,TDATE,DTR,RTD,
  6.      2TPREC,PI,P,Q,V,U,A1,A2,A3,A4,A5,A6
  7.       COMMON /ASTRO/OELMT,TDATE,DTR,RTD,PI,CLAT,CRAD
  8.       TPREC = (TDATE-2415020.0D0)/36525.0D0
  9. C    ELEMENTS FOR EARTH
  10.       OELMT(1,2) = TDATE
  11.       OELMT(2,2) = ((((1.4D-8*TPREC)-9.33D-6)*TPREC)+1.30855D-2)*TPREC-
  12.      26.54D-3
  13.       OELMT(3,2) = ((((-2.8D-8*TPREC)+7.94D-6)*TPREC)-2.415735D-1)*
  14.      2TPREC+174.52817D0
  15.       OELMT(4,2) = ((((3.333D-6*TPREC)+1.361D-4)*TPREC)+5.647073D-1)*
  16.      2TPREC+287.390758D0
  17.       OELMT(7,2) = ((-1.26D-7*TPREC)-4.18D-5)*TPREC+1.675104D-2
  18.       OELMT(8,2) = ((-3.3D-6*TPREC)-1.5D-4)*TPREC**2
  19.       OELMT(8,2) = OELMT(8,2)+DMOD((TPREC*3.599904975D4+3.5847583D2)
  20.      2,360.0D0)
  21.       A1 = DTR*(153.23+TPREC*2.25187541D4)
  22.       A2 = DTR*(216.57+TPREC*4.50375082D4)
  23.       A3 = DTR*(312.69+TPREC*3.29643577D4)
  24.       A4 = DTR*(350.74+TPREC*(4.452671142D5-(TPREC*1.44D-3)))
  25.       A5 = DTR*(231.19+TPREC*2.02D1)
  26.       A6 = DTR*(353.4+TPREC*6.59287155D4)
  27.       CLAT(2) = DTR*(1.34D-3*COS(A1)+1.54D-3*COS(A2)*2.0D-3*COS(A3)+
  28.      21.79D-3*SIN(A4)+1.78D-3*SIN(A5))
  29.       CRAD = 5.43D-6*SIN(A1)+1.575D-5*SIN(A2)+1.627D-5*SIN(A3)+
  30.      23.076D-5*COS(A4)+9.27D-6*SIN(A6)
  31. C    ELEMENTS FOR VENUS
  32.       OELMT(1,3) = TDATE
  33.       OELMT(2,3) = ((-3.251D-5*TPREC)-8.226D-4)*TPREC+3.394552D0
  34.       OELMT(3,3) = ((-1.4039D-4*TPREC)-2.777139D-1)*TPREC+76.368593D0
  35.       OELMT(4,3) = ((-1.14435D-3*TPREC)+2.893249D-1)*TPREC+54.493527D0
  36.       OELMT(7,3) = ((9.1D-8*TPREC)-4.774D-5)*TPREC+6.82069D-3
  37.       OELMT(8,3) = ((3.097D-4*TPREC)+5.851921191D4)*TPREC+342.767053D0
  38.      2-(OELMT(4,3)+OELMT(3,3))
  39. C    ELEMENTS FOR MARS
  40.       OELMT(1,4) = TDATE
  41.       OELMT(2,4) = ((-2.305D-5*TPREC)-8.1839D-3)*TPREC+1.854113D0
  42.       OELMT(3,4) = ((-6.4435D-4*TPREC)-2.940497D-1)*TPREC+49.319212D0
  43.       OELMT(4,4) = ((4.6647D-4*TPREC)+7.385934D-1)*TPREC+285.597172
  44.       OELMT(7,4) = ((-7.7D-8*TPREC)+9.2064D-5)*TPREC+9.33129D-2
  45.       OELMT(8,4) = ((3.107D-4*TPREC)+1.914169551D4)*TPREC+293.
  46.      2737334D0-(OELMT(4,4)+OELMT(3,4))
  47. C    ELEMENTS FOR JUPITER
  48.       OELMT(1,5) = TDATE
  49.       OELMT(2,5) = ((2.952D-5*TPREC)-2.2192D-3)*TPREC+1.307028D0
  50.       OELMT(3,5) = ((9.5857D-4*TPREC)+1.661852D-1)*TPREC+99.865881D0
  51.       OELMT(4,5) = ((-2.1041D-4*TPREC)+4.7591D-2)*TPREC+273.553214D0
  52.       OELMT(7,5) = ((-4.676D-7*TPREC)+1.6418D-4)*TPREC+4.833475D-2
  53.       OELMT(8,5) = ((3.347D-4*TPREC)+3.036301986D3)*TPREC+238.049257D0
  54.      2-(OELMT(4,5)+OELMT(3,5))
  55.       V = TPREC/5.0D0+1.0D-1
  56.       P = 237.47555+TPREC*3034.9061
  57.       Q = 265.9165+TPREC*1222.1139
  58.       U = DTR*(5.0*Q-2.0*P)
  59.       CLAT(5) = DTR*(.331364-V*1.0281D-2-V*V*4.692D-3)*SIN(U)
  60. C    ELEMENTS FOR SATURN
  61.       OELMT(1,6) = TDATE
  62.       OELMT(2,6) = ((-5.022D-5*TPREC)+2.419D-3)*TPREC+2.489374D0
  63.       OELMT(3,6) = ((-1.8862D-4*TPREC)-2.597237D-1)*TPREC+113.356715D0
  64.       OELMT(4,6) = ((7.0612D-4*TPREC)+8.218494D-1)*TPREC+338.439665D0
  65.       OELMT(7,6) = ((-4.676D-7*TPREC)+1.6418D-4)*TPREC+4.833475D-2
  66.       OELMT(8,6) = ((3.245D-4*TPREC)+1.22350988D3)*TPREC+2.66564377
  67.      2D2-(OELMT(4,6)+OELMT(3,6))
  68.       V = TPREC/5.0D0+1.0D-1
  69.       P = 237.47555+TPREC*3034.9061
  70.       Q = 265.9165+TPREC*1222.1139
  71.       U = DTR*(5.0*Q-2.0*P)
  72.       CLAT(6) = DTR*(-.814181D0+V*1.815D-2+V*V*1.6714D-2)*DSIN(U)
  73. C    ELEMENTS FOR URANUS
  74.       OELMT(1,7) = TDATE
  75.       OELMT(2,7) = ((-2.2D-7*TPREC)-1.7599D-3)*TPREC+7.73723D-1
  76.       OELMT(3,7) = ((4.2988D-4*TPREC)+5.57505D-2)*TPREC+73.700227D0
  77.       OELMT(4,7) = ((-5.0125D-4*TPREC)+3.2554D-2)*TPREC+98.546561D0
  78.       OELMT(7,7) = ((7.7D-8*TPREC)-2.658D-5)*TPREC+4.63444D-2
  79.       OELMT(8,7) = ((3.16D-4*TPREC)+4.29863546D2)*TPREC+2.441974D2-
  80.      2(OELMT(4,7)+OELMT(3,7))
  81. C    ELEMENTS FOR NEPTUNE
  82.       OELMT(1,8) = TDATE
  83.       OELMT(2,8) = ((-2.27D-6*TPREC)-1.5D-5)*TPREC+1.774485D0
  84.       OELMT(3,8) = ((4.435D-5*TPREC)-8.3952D-3)*TPREC+131.234637D0
  85.       OELMT(4,8) = ((3.842D-5*TPREC)+3.67891D-2)*TPREC+276.190852D0
  86.       OELMT(7,8) = ((-2.0D-9*TPREC)+6.33D-6)*TPREC+8.99704D-3
  87.       OELMT(8,8) = ((3.205D-4*TPREC)+2.19555914D2)*TPREC+84.457994D0
  88.      2-(OELMT(4,8)+OELMT(3,8))
  89.       RETURN
  90.       END
  91.